Medium Modification of 7-jets in High-energy Heavy-ion Collisions 
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Medium modification of 7-tagged jets in high-energy heavy-ion collisions is investigated within 
a Linearized Boltzmann Transport model for jet propagation that includes both elastic parton 
scattering and induced gluon emission. Inclusion of recoiled medium partons in the reconstruction 
of partonic jets is found to significantly reduce the net jet energy loss. Experimental data on 7-jet 
asymmetry and survival rate in Pb + Pb collisions at yfs = 2.76 TeV can be reproduced. Medium 
modifications of reconstructed jet fragmentation function, transverse profile and energy flow outside 
the jet-cone are found to be sizable especially for 7-tagged jets with small values of x = p^/Pr- 
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Jet quenching [ij, caused by jet- medium interaction 
and manifested in the suppression of large P t hadron 
spectra, dihadron and 7-hadron correlation, has been 
used successfully to probe properties of quark-gluon 
plasma in high-energy heavy-ion collisions. The study 
has been extended to reconstructed jets which are less 
sensitive to the non-perturbative process of hadroniza- 
tion. Large dijet asymmetry in central Pb + Pb colli- 
sions at y/s — 2.76 TeV is indeed observed at the Large 
Hadron Collider (LHC) @]. Further studies of experi- 
mental data within models of jet propagation 0-01 can 
yield information on the strength of jet-medium interac- 
tion and transport properties of the quark-gluon plasma. 

A jet in the partonic picture contains many shower 
partons from the final-state radiation of the leading 
parton. These shower partons will interact with the 
thermal medium and lead to parton energy loss, trans- 
verse momentum broadening and modification of the jet 
shape. Transport of the recoiled thermal partons from 
jet-medium interaction, on the other hand, can also lead 
to jet- induced medium excitation Q. Jets in heavy- ion 
experiments are reconstructed with a jet-finding algo- 
rithm [9| and consist of collimated clusters of hadrons 
within a jet-cone y/ (cj> — 4>c) 2 + (v — Vc) 2 < R } where 77 
and <f> are hadrons' pseudo-rapidity and azimuthal angle, 
respectively, and the subscript C denotes the center of 
the jet-cone. Energy within the jet-cone should contain 
hadrons from both shower and recoiled thermal partons. 
It is therefore important to analyze experimental data 
taking into account of both medium- modified shower par- 
ton distributions and jet-induced medium excitation. 

Back-to-back 7-jets are considered "golden channels" 
for the study of jet quenching since they are free of trigger 
bias [Io| as compared to dijets. Though 7-jet asymmetry 
as measured in Pb + Pb collisions at LHC [ill Il3 can be 



ification of 7-tagged jets in high-energy heavy-ion colli- 
sions within a Linearized Boltzmann Transport (LBT) 
model [l5[ which tracks the transport of both shower and 
recoiled medium partons. We will focus particularly on 
the effect of jet-induced medium excitation on jet energy 
loss, modification of jet profiles and energy flow outside 
the jet-cone. We further study the dependence of jet 
profile on jet propagation length by selecting jets with 
different 7-jet asymmetry. 

Wit hin the LBT model, propagation of jet shower par- 
tons and medium excitation is simulated according to a 
linearized Boltzmann equation, 
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dfl{pi) = - J dp 2 dp 3 dp 4 (f 1 f2 - /3/4)|Ml 2 ->34| 
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explained by parton energy loss |13l Il4| , no attention has 
been paid to the effect of jet-induced medium excitation. 
In this Letter, we will report a first study of medium mod- 



where d Pi = d 3 p l /[2S I (2 7 r) 3 ], /< = l/( e f"/ T ± 1) (i = 
2, 4) are parton phase-space distributions in a thermal 
medium with local temperature T and fluid velocity u = 
(l,v)/VT=tf, fi = (2n) 3 5 3 (p- Pl )6 3 (x-x t - v\t) (i = 
1,3) are the parton phase-space densities before and after 
scattering. We assume small angle approximation for 
the elastic scattering amplitude |Mi2^34| 2 = C<? 4 (s 2 + 
u 2 )/(—i + /sjj) 2 with Debye screening mass hd, where 
s, t, and u are Mandelstam variables, and C=l (9/4) is 
the color factor for quark-gluon (gluon-gluon) scattering. 
The strong coupling constant a s = <? 2 /47r is fixed and 
will be determined via comparison to experimental data. 

Partons are assumed to propagate along classical tra- 
jectories between two adjacent collisions. The probabil- 
ity of scattering is determined in each time step At by 
P a = 1 - exp[- Y,j(^ x j ' u ) Eb VabPbixj)}, where a ab is 
the parton scattering cross section, the sum over time 
steps starts from the last scattering point, and p b is the 
local medium parton density. Both shower (^3) and re- 
coiled medium partons (^4) after each scattering are fol- 
lowed by further scatterings in the medium. To account 
for the back-reaction in the Boltzmann transport, initial 



2 



thermal partons (jp2), denoted as "negative" partons, are 
transported according to the Boltzmann equation. Their 
energies and momenta will be subtracted from all final 
observables. These "negative" partons are considered as 
part of the recoiled partons that are responsible for jet- 
induced medium excitations (IBl ]. 

In this study, the LBT is extended to include induced 
radiation accompanying each elastic scattering according 
to the high- twist approach flil ]. 
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where z and kj_ are the energy fraction and transverse 
momentum of the radiated gluon, p^ — p^/po, P(z) = 
[l + (l-z) 2 ]/z the splitting function, 77 = 2Ez{l-z)/k\ 
the gluon formation time, and q a = ^2 h Pb J dtq\da a bl dt 
the jet transport parameter. The Debye screening mass 
fin is used as an infrared cut-off for the gluon's energy. 
Multiple gluon emission induced by a single scattering is 
included according to a Poisson distribution. All radiated 
gluons are assumed to be on-shell and their 4-momcnta 
are successively determined from Eq. @. The first ra- 
diated gluon is assumed to be associated with the cur- 
rent scattering with given transverse momentum transfer. 
The 4-momenta of the final beam and recoiled parton are 
determined according to a 2 — > 3 process. The next gluon 
emission is carried out by setting the final beam par- 
ton of the previous gluon emission off-shell with its large 
component of the light-cone momentum unchanged. The 
transverse momentum of the previous emitted gluon will 
be reset by its on-shell condition and 4-momentum con- 
servation of the current gluon emission. Such procedure 
can be repeated for the rest of N g gluon emissions. 

For initial configurations of 7-jets, we use HIJING [13] 
for p + p collisions at yfs = 2.76 TeV with a trigger on 
the transverse momentum transfer qr > 30 GeV in the 
CM frame of two colliding partons. Further selections are 
made for events with p7p > 60 GeV. Jet shower partons, 
including those from both initial and final-state radia- 
tion, are transported through a thermal medium within 
the LBT model. The final partons, including "negative" 
partons, are used for jet reconstruction using a modi- 
fied version of the anti-fc t algorithm in FAST JET in 
which energies and momenta of "negative" partons are 
subtracted from the final jet observables. Energies of 
jets reconstructed from final hadrons and partons differ 
only about 1 GeV in p + p collisions from HIJING. 

To illustrate the influence of jet-induced medium exci- 
tation on the final jet energy, we consider first the prop- 
agation of 7-tagged jets in a uniform and static gluonic 
medium at temperature T = 300 MeV. We set a s = 0.4 
and n 2 D = 1 GeV 2 for simplicity and use a jet-cone size 
R = 0.4. Shown in Fig. QJa) is the net jet energy loss as 
a function of the propagation time. If we include only jet 
shower partons (dotted) in the jet reconstruction, the en- 
ergy loss is considerably larger than when radiated gluons 
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FIG. 1: (a) Energy loss as a function of time; (b) azimuthal 
distribution relative to the 7-direction, (c) reconstructed jet 
fragmentation function (parton) at different times, (d) initial 
(dotted) and final jet transverse profile with (solid) and with- 
out (dashed) induced radiation for 7-tagged jets in a uniform 
gluonic medium. See text for more detailed descriptions. 



(dashed) or all (shower + radiated + recoiled ) partons 
(solid) arc included. Therefore, both radiated gluons and 
recoiled partons from the elastic scattering enter the jet- 
cone and become part of the reconstructed jets. The net 
jet energy loss shows a clear quadratic time-dependence 
during the early times. The reconstructed jets also show 
significant deflection with respect to the 7-direction at 
later times due to jet-medium interaction as shown in 
Fig.QJb). Shown in Fig. [He) is the longitudinal momen- 
tum fraction (z = pL/E }ct ) distribution, which we will 
refer to as reconstructed jet fragmentation function, of 
partons within the jet-cone at different times. There is 
very little change at large and intermediate z, but huge 
enhancement at low z due to contributions from recoiled 
partons and radiated gluons. In Fig. QJd) we show the 
jet transverse profile, 
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of the initial jet (dotted) and at a later time t = 8 fm/c 
(solid and dashed), where pt^i,?^) is the summed pr 
of all partons in the annulus between radius r\ and 
inside the jet-cone. We also extend the above transverse 
energy profile to outside the jet-cone between R < r < 
it/2. Without induced gluon radiation, elastic scatterings 
have little effect on the jet transverse profile within the 
jet-cone while recoiled partons contribute to transverse 
energy flow outside the cone. Induced gluon radiation, 
on the other hand, significantly broadens the transverse 
energy profile both within and outside the jet-cone. 



3 



t (fm/c) 
10 



0.2 



I 



( b ) 0-10% PbPb 
Vs=2.76 TeV 
• CMS 




10 



- 10 



10 



f 

-a 
-l 



10 



10 



FIG. 2: The same as Fig. Q] except in central (0-10%) Pb + 
Pb collisions at yfs — 2.76 TeV. The inserts in (c) and (d) 
are ratios of initial and final reconstructed jet fragmentation 
functions (parton) and transverse profiles, respectively. See 
text for detailed descriptions. 



To study the effect of jet energy loss and medium mod- 
ification of 7-tagged jets in heavy-ion collisions within 
LBT model, we use the space-time profile of tempera- 
ture and fluid velocity in the quark-gluon phase from 
(3+l)D ideal hydrodynamical simulations 18j of Pb + 
Pb collisions at the LHC. The initial 7-jet production 
from HIJING is distributed according to the overlap func- 
tion of two colliding nuclei with a Wood-Saxon nuclear 
geometry. Whenever comparisons are made to the exper- 
imental data, we apply the same kinematic cuts to LBT 
results. 



For CMS data 11 
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> 60 GeV, 
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and for ATLAS data [3, 60 < V t < 90 GeV 



rf\ < 1.44, 
0T| > 7tt/8; 
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< 1.3, 



rf? > 25 GeV, 



jet 



< 2.1, and A0 > 7tt/8. In LBT sim- 



ulations, we use a Debye screening mass y? D = Aira s T 2 . 

Shown in Fig. Eta) is the averaged jet energy loss as a 
function of time in the most 10% central Pb + Pb colli- 
sions at s/s = 2.76 TeV for a jet cone-size R = 0.3 from 
LBT simulations with a s = 0.2. Similar to the results in 
a uniform medium in Fig. QJa), inclusion of recoiled par- 
tons (solid) significantly reduces the net jet energy loss 
as compared to the case when only shower partons (dot- 
ted) and radiated gluons (dashed) are included in the jet 
reconstruction. Because of rapid expansion and cooling 
of the thermal medium, jet energy loss saturates approx- 
imately 10 fm/c after an initial linear rise. For the same 
reason the jet azimuthal distribution as shown Fig. [2jb) 
remains almost unchanged as compared to that in p + 
p collisions, in agreement with CMS data [ll|. Modifi- 
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FIG. 3: Distribution of 7-jet asymmetry x = PjIVPt in Pb 
+ Pb collisions with four centralities at y/s — 2.76 TeV from 
LBT with a B = 0.2 as compared to CMS data [ll[. 



Fig- [U[c) is also very small at large and intermediate z 
region but sees big enhancement at small z due mostly 
to contributions from radiated gluons. The jet transverse 
profile shown in Fig.[5Jd) is enhanced by 40% around the 
edge both inside and outside the jet-cone. The jet trans- 
verse profile outside the jet-cone are also enhanced due 
to dissipation of the lost jet energy that are carried by re- 
coiled partons and radiated gluons. Similar features have 
been observed in single jets in central Pb + Pb collisions 
at LHC 0. 

Following experiments at LHC [H, 12|, one can study 
the effect of jet energy loss on the 7-jet asymmetry distri- 
bution, dN/dx with x = p^/Pr- Shown in Fig. [3] are 7- 
jet asymmetry distributions from LBT simulations (his- 
togram) as compared to CMS data [ll[ in Pb + Pb col- 
lisions with four different centralities. Because of initial- 
state radiation, 7-jets are produced with large momen- 
tum asymmetry in both p + p and peripheral heavy-ion 
collisions where jet energy loss is small. The asymmetry 
distributions are, however, skewed to smaller values of x 
in central heavy-ion collisions due to jet energy loss as 
seen in the LBT results which can fit the experimental 
data quite well with a fixed value of a s = 0.2. 

We can further quantify the 7-jet asymmetry in heavy- 
ion collisions by the averaged asymmetry or ratio of the 
jet and photon transverse momenta (x) = (p^ /Pj^j an d 
the jet survival rate or fraction of 7-tagged jets i?j 7 with 
p^ 4 > 30 GeV (CMS cut) or x = p^/Pr > °- 42 (AT- 
LAS cut). Shown in Fig. [4] are LBT results (lines) on 
the averaged 7-jet asymmetry (x) and jet survival rate 
Rjj as functions of the number of participant nucleons 



cation of the reconstructed jet fragmentation function in 



are compared to CMS (solid circle and square) [llj and 
ATLAS (open circle) data [l~2| ■ Note that kinematic cuts 
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FIG. 4: Averaged 7-jet asymmetry (x) = {p 1 ^ IPtJ (left) 
and jet survival rate Rj 1 (right) as functions of the number 
of participant nucleons in Pb+Pb at y^s = 2.76 TeV from 
LBT as compared to experimental data [Til, ■ Values of 
a s = 0.15 - 0.23 (dashed) and 0.2 - 0.27 (solid) are used for 
LBT calculations with CMS and ATLAS cuts, respectively. 
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the LBT model that can reproduce well recent exper- 
imental data on 7-jet asymmetry and survival rate in 
Pb + Pb collisions at LHC. Inclusion of recoiled medium 
partons in the jet reconstruction is found to lead to signif- 
icant reduction in net jet energy loss therefore should be 
taken into account in model studies. No jet deflection are 
observed and medium modification of the reconstructed 
jet fragmentation function and transverse profile are also 
quite modest because of the short averaged path length 
in an expanding system. By selecting events with small 
values of x = p}^ /pX, one can focus on jets with large 
energy loss and therefore will see much stronger medium 
modification. Such approach can help to achieve more 
detailed jet tomography of the quark-gluon plasma in 
high-energy heavy-ion collisions. 
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FIG. 5: Medium modification of the reconstructed jet frag- 
mentation function (parton) (left) and transverse profile 
(right) with momentum asymmetry x < 0.5 (solid) and 
x > 1.0 (dashed) for 7-tagged jets in central (0-10%) Pb + 
Pb collisions at Js = 2.76 TeV. 



in ATLAS data, which are also truncated with x < 0.42, 
are somewhat different from that in CMS data. In LBT 
calculations, a s =0. 15-0.23 are used with the CMS cuts 
(dashed) while a s =0.2-0.27 for the ATLAS cuts (solid). 
One can see that the averaged momentum asymmetry 
(x) has a very weak centrality dependence and is not very 
sensitive to the value of a s . The jet survival rate i?j 7 has, 
however, a stronger dependence on the centrality and the 
value of a s or strength of jet-medium interaction. 

Since smaller values of x = /Pt imply larger aver- 
age jet energy loss, one should expect stronger medium 
modification of 7-tagged jets. This is indeed the case as 
seen in Fig. [5j where we show medium modification of 
the reconstructed jet fragmentation function for partons 
(left) and transverse profile (right) for 7-tagged jets with 
different values of x. The reconstructed jet fragmenta- 
tion function at low z and transverse profile toward the 
edge of the jet-cone are both significantly enhanced for 
events with x < 0.5. One should expect to see larger en- 
hancement for smaller values of x which can be achieved 
by increasing the energy of the triggered 7. 

In summary, we have studied medium modification of 
7-tagged jets in high-energy heavy-ion collisions within 
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